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SUMMARY 

This paper presents a computer-aided design procedure for generating bevel 
gears. The development is based on examining a perfectly plastic, cone-shaped 
gear blank rolling over a cutting tooth on a plane crown rack. The resulting 
impression on the plastic gear blank is the envelope of the cutting tooth. 

This impression and envelope thus form a conjugate tooth surface. Equations 
are presented for the locus of points on the tooth surface. The same proce- 
dures are then extended to simulate the generation of a spiral bevel gear. 

The corresponding governing equations are presented. 


INTRODUCTION 

Bevel gears are the principle means of motion and power transmission 
between intersecting shafts. Their use is widespread. The geometrical cha- 
racteristics of bevel gears have long been documented by the American Gear 
Manufacturers' Association (AGMA, 1964) and others (Dudley, 1962; Dyson, 1969; 
Bonsignore, 1976; Litvin et al., 1975, 1982, 1983; Baxter, 1966; Huston and 
Coy, 1 981 , 1 982a, b) . Recent advances in computer-aided design present opportu- 
nities for a new look at the geometry of these gears. These computer-based 
procedures also provide a means for optimizing the geometry. 

In a previous paper (Chang, Huston, and Coy, 1984) we have demonstrated 
the feasibility of this procedure to determine a spur gear tooth profile. The 
basic idea was to use the envelope of a family of curves to develop an invo- 
lute spur gear tooth profile. The study demonstrated that the envelope of an 
inclined straight-line segment on a rolling gear blank is an involute of a cir- 
cle. The inclined line segment in turn represented the rack tooth of a hob 
cutter. Such a procedure provides a means for analytically and numerically 
determining the tooth profile, given a cutter profile. 


*Work funded under NASA Grant NAG3-188. 



In the current paper we extend these ideas to the generation of both 
straight and spiral bevel gears. The paper itself is divided into four sec- 
tions with the section following the Symbols providing preliminary ideas use- 
ful in this sequel. The next two sections describe the formulation of string 
and spiral bevel gear tooth surfaces. This is followed by a discussion of 
appl i cations . 

SYMBOLS 

The following symbols are used in the section for development of a 
straight bevel gear tooth. 

R m Mean radius of crown gear in pitch plane 

[R] Coordinate transformation matrix; from S to S 

[R-|] Coordinate transformation matrix; from S] to S 

C R£ ] Coordinate transformation matrix; from S2 to S] 

[R3] Coordinate transformation matrix; from S to S2 

r Position vector in S 

r g Position vector in S 

r i j Components of [R] 

S(X,Y,Z) Cartesian coordinate system stationary with ground 

5 (X,Y,Z) Cartesian system stationary with generated gear 

S-|(X-| ,Yq ,Z-[) Intermediate Cartesian coordinate system 
S2<X2,Y2,Z2) Intermediate Cartesian coordinate system 

Xp Projection of X axis on crown gear pitch plane 

a Angular position of gear referenced to system S 

aQ Initial angular position of gear 

Y Pitch angle of pinion 

6 Angle of rotation of generated gear 

4> Complement of pressure angle of cutter 

Y Cutter's orientation angle on pitch plane 

The following symbols are used in the section for development of a spiral 
bevel gear tooth. 

B Derivative of T with respect to $2 


2 



b i 

d U 

CE] 

e u 

H 

h 

Ng 

N P 

[R c i] 

[R 1g ] 

CRgp] 

[R p2 ] 


n 

^ 2 

Sc(X c ,Yc,Z c ) 
S i < X i , Y ! , Z -| ) 
S 2 (X2,Y2,Z 2 ) 

Sg(Xg , Yg , Zg) 
Sp(Xp ,Yp , Zp) 


Components of B 
9e i j / 8<j> 2 

Coordinate transformation matrix; from S 2 
Components of CE] 

Horizontal machine setting of cutter 
Addendum of cutter tooth (pitch to tip) 
Number of teeth in crown gear 
Number of teeth in generated gear 
Coordinate transformation matrix; from S] 
Coordinate transformation matrix; from Sg 
Coordinate transformation matrix; from Sp 
Coordinate transformation matrix; from S 2 
Mean radius of cutter in pitch plane 
Position vector in S c 
Position vector in Sg 
Position vector in Sp 
Position vector in Si 
Position vector in S 2 

Cartesian coordinate system stationary with 


to S c 


to S c 
to Si 

tO Sg 

to Sp 


cutter with origin 


at 0 C 

Cartesian coordinate system stationary with gear with origin 
at Oi 

Cartesian coordinate system stationary with pinion with 
origin at 0 2 

Cartesian coordinate system stationary with ground with 
origin at Og 

Cartesian coordinate system stationary with ground with origin 
at Op 
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Coordinate translation transformation matrix; from S 2 to S c 
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Coordinate translation transformation matrix; from S c to Si 

Coordinate translation transformation matrix; from Sp to Sg 
Components of T 

Vertical machine setting of cutter 
Angular speed ratio 
Z c intercept of cutter surface 
Surface coordinate of cutter 
Root angle of pinion 
Pitch angle of pinion 

Angle of rotation of system S c reference to S] 

Angle of rotation of gear 
Angle of rotation of gear blank 
Initial angular position of gear 
Complement of pressure angle of cutter edge 

Rotation rate of gear 
Rotation rate of pinion 


ENVELOPE OF A FAMILY OF SURFACES 

Consider a family of surfaces represented by an equation of the form 
F< x i ,X 2 ,X 3 ,t) = 0 with a parameter t. Let S be a surface of the family and 
let it be intersected by neighboring surfaces S'. If S and S' correspond 
to the values t and t + At, the curve is represented by the simultaneous 
equations 


F ( x i ,X2,X3,t) = 0 \ 

F( x i ,X2,X3,t + At) = 0 J 


It may also be represented by the equations 


F ( x i ,X2,X3,t) = 0 

F(x,,X2,x t + At) - F(x 1 .Xg.x^.t) 


At 


= 0 


( 2 ) 
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The surface [F(X] ,X 2 ,X 3 ,t + At) - F(x] ,X 2 ,X 3 ,t)]/At = 0 goes through the curve 
common to the two surfaces F(x] ,X 2 ,X 3 ,t) = 0 and F(xi,X 2 ,X 3 > t + At) = 0. 

When S' approaches S as a limit (i.e., when At approaches zero) the inter- 
section curve will approach a limiting curve. This curve is given by 


F( x i ,X2,X3,t) = 0 

^.x^.t) = 0 


(3) 


Equation (3), when t is fixed, represents a curve on the 
ily. The same equation, with t variable, will represent 
and will generate a surface. This surface is the envelope 
ily. The result of eliminating t is the equation of the 
1935). 


surface of the fam- 
a fami Ty of curves 
of the given fam- 
envelope (Graustein, 


In the two-dimensional case the envelope of a family of curves is a curve tan- 
gent to the given family. For example, the envelope given by an inclined 
straight-line segment on a rolling gear blank, is found to be an involute of a 
circle (Chang, Huston, and Coy, 1984). 

In the following sections these procedures are generalized to simulate the 
straight bevel and spiral bevel gear tooth surface generation process. The 
family of surfaces created by the cutter generates an envelope in the gear 
blank. The envelope in the gear blank then forms a conjugate tooth surface. 


DEVELOPMENT OF A STRAIGHT BEVEL GEAR TOOTH SURFACE 

The cutter used for the surface generation of a straight bevel gear tooth 
is called the basic crown rack. Figure 1 depicts the machining model of the 

straight bevel gear tooth generation process. The basic crown rack R is con- 

sidered to be fixed on the imaginary crown gear pitch plane C as a rack step, 
as depicted in Fig. 2. The generated gear blank G is a cone with the vertex 
at the machine center 0. The gear blank is allowed to roll on the crown gear 
pitch plane C. The crown gear pitch plane is an imaginary fixed plane. When 
the gear blank rolls over the crown rack step, the envelope of the basic crown 

rack forms the tooth surface of the straight bevel gear. 

The coordinate system used to describe^the crown gear is S(X,Y,Z,) with 
origin at 0. The coordinate system S(X,Y,Z) is fixed on the gear blank with 
origin at 0. The gear blank is allowed to rotate through an angle 9 about 
X. The pitch angle of the gear blank is denoted as y. The initial position 
of the gear blank is defined by angle aQ, which is the angle between the 
axis X and^the projection of the X axis on the crown gear pitch plane, 
denoted by Xp. During the gear blank rolling motion the angular position of 
the gear blank is defined by angle a, where a is measured between Xp and 
X. From this geometry we obtain the relation between 0 and a as 


a = aQ - 9 tan y (4) 

The coordinate system S may be obtained by coordinate transformation of 
the system S in the following steps: 


5 



Step 

Axis of 
rotation 

Angle 

turned 

Coordinate 

system 

Rotation 

matrix 

0 

— 

— 

At XYZ initially 

— 

1 

Z 

a 

XlYiZi 

R 1 

2 

Yl 

-Y 

x 2 y 2 z 2 

r 2 

3 

*2 

e 

XYZ 

r 3 


The indicated rotation matrices are 


cos a -sin a 


CR, ] = 

sin 

a 
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a 
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Y 
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0 

sin 
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(5) 


( 6 ) 


(7) 


The posi tion^vectors r with the components (x,y,z) and^ r g with the 
components (x,y,z) locating a typical point in S and S are related by the 
expression 

r = [R]] [R 2 ] C R 3 ] r g = [R]r g (8) 

From Eqs. (5) to (7) the elements of [R], r^(i,j = 1,3), are 


\ 


ri 1 * COS a COS Y 

r 12 = -cos a sin y sin 6 - sin a cos 9 

r ] 3 = -cos a sin y cos 6 + sin a sin 0 

r 2 i = sin a cos y 

r 22 = -sin a sin y sin 9 + cos a cos 0 

r 23 = -sin a Y cos 7 * 9 - cos a sin 9 

r 3 1 = sin y 

r3 2 = sin 0 cos y 

r33 = cos 0 cos y 


( 9 ) 




i 
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EQUATION FOR CROHN RACK AND STRAIGHT BEVEL GEAR TOOTH SURFACE 

The basic crown rack may be viewed as a straight cutting blade reciprocat- 
ing in the radial direction. The tooth surface of the rack thus forms an 
inclined plane passing through the crown gear center. The pitch plane of the 
cutter surface is shown in Fig. 2. The normal plane view in Fig. 3 shows the 
rack profile. The equation of the rack surface may be expressed in terms of 
system S as 


x tan i|> - y - z cot 4> = 0 (10) 

where y is the angle between the cutter face and the XZ plane, and <j> is 

the complement of the pressure angle of the cutter as shown in Fig. 2. 

The cutter's surface may be expressed in terms of the gear blank's system S 

by substituting Eqs. (8) and (9) into Eq. (10). This leads to 

x[(tan y cos a - sin <x)cos y - cot <j> sin y] 

+ y[sin y sin 9(sin a - tan ^ cos a) - cos 6(cos a + tan y sin a) 

- cot <j> cos y sin 0] 

+ zCsin y cos 0(sin a - tan y cos a) - sin 0(cos a + tan y sin a) 

- cot cos y cos 0] = 0 = F(x,y,z,0) (ID 

Following this procedure, to determine the envelope of the cutter surface 
on the generated gear blank, we evaluate the partial derivative of Eq. (11) 
with respect to parameter 0, the rotation angle of the gear blank. This 
produces the relation 

x[sin y(tan <J> sin a + cos a)] 

+ ytsin 0(1- sin y tan yXcos a + tan y sin a) 

+ cos 0(sin y - tan yXsin a - tan cos a) - cot <|> cos y cos 0] 

+ z[cos 9(1 - sin y tan yXcos a + tan y sin a) 

- sin 0(sin y - tan yXsin a - tan cos a) + cot <f> cos y sin 0] = 0 (12) 

Let x be an independent variable. Then, from Eqs. (11) and (12), y 

and z may be evaluated in terms of x. Observe that since the coefficients 
are functions of 9, the tooth surface has the parametric form 

/-V /-V \ 

X = X \ 

y * y(x ,0) / (13) 

z = z(x ,9) / 
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where x and 0 are the surface coordinates. Equation (13) represents the 
envelope of the rack relative to the gear blank. This represents the tooth 
surface of the straight bevel gear. 


DEVELOPMENT OF A SPIRAL BEVEL GEAR TOOTH SURFACE 

A spiral bevel gear tooth surface is developed in a similar manner. Fig- 
ure 4 depicts a circular cutter generating a spiral bevel gear. The cutter is 
mounted on the cradel of a generating machine. When the cutter rotates about 
its own axis, it forms a surface that simulates a crown gear. As the cradle, 
and hence the cutter, rotates about Zg at the rate cog and the gear blank 
rotates about Z p at the rate w p , the cutter will generate a spiral bevel 
gear. The cutting speed is independent of cog and co p . It is not related to 
the kinematics of tooth generation. The relation between o>g and co p is 
simply 


Wp/Wg = Ng/N p 


(14) 


where Ng and N p are the numbers of teeth in the crown gear and generated 
gear, respectively. 

The coordinate system used to describe the crown gear is Si<Xi,Y],Zi> 
with origin at 0]. The crown gear G with frame S] fixed in G rotates 
through an angle <J>] about Zg with respect to a global coordinate system 
Sg(Xy, Y g,Zg) with the origin at 0g as_in Fig. 4. The position vectors r] , 
with the components (x],y],Z]), and rg, with the components (Xg.yg.Zg), locat- 
ing a point in Si and Sg are related by the expression (see y Fig. 5) 


H - ^lg^ ^g 


(15) 


where [Rig] is an orthogonal transformation matrix given by 


[R,g] 


cos 4> 1 sin 4> 1 0 

-sin <j>j cos <j> 1 0 

0 0 1 


(16) 


Let a coordinate system S C (X C ,Y C ,Z C ) be fixed on the cutter with origin 
at 0 c (H,V,0). Let H and V be the horizontal and vertical machine settings 
(see Fig. 5). The cutter rotates through an angle 0 about axis Z c . Posi- 
tion vectors r c , with the components (x c ,x c ,z c ) locating a point relative to 
S c , and n are related by the expression 


r c ■ tR cl ] r l + 


[R d ] T lc 


(17) 
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where CR c i 3 and Ti c are 


cos 6 sin e 0 

CR cl ] = -sin e cos 0 0 (18) 

0 0 1 

and 

H 

T ,c ' V <’ 9) 

0 

Let the coordinate system S 2 (X 2 ,Y 2 ,Z 2 ) be fixed in the gear blank. Let 
S 2 rotate through an angle <j >2 about Zp with respect to a second global 
coordinate system Sp(X p ,Yp,Zp) (see Figs. 4, 6, and 7). The position vec- 
tors r p , with the components (Xp,y p ,Zp), and r 2 , with the components 
(x 2 .y 2 .z 2 >> locating a point in S p and S 2 are related by the expression 

r p = [R p2 1 ?2 <20) 

where the transformation matrix CR p 2 3 from S 2 to S p is given by 

cos <t >2 -sin 4 > 2 0 

CRp 2 3 = sin <J> 2 cos 4> 2 0 (21 ) 

0 0 1 

The two global coordinate systems Sg and S p are related by the root 
angle y of the generated gear and the addendum of the cutter tooth h (see 
Figs. 4, 6, and 8.) Hence, rg and r p are related by the expression 

fg = tRgp]r p + Tg p (22) 

where 

1 0 0 

C Rgp 3 = 0 sin y -cos y (23) 

0 cos y sin y 

and 


9 


( 24 ) 
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During the cutting process the simulated crown gear rotates in such a way 
that the motion is conjugate with the generated gear blank. In Fig. 6 the 
pitch element O 2 P is an instantaneous axis for these "conjugate gears." 

Hence, their angular velocity components on the pitch element are equal. That 
is, 


<i>2 sjn Y 0 = <M COS(y 0 - Y> 


(25) 


where yo is the pitch angle of the generated gear. Integrating Eq. (25) 
with respect to time then leads to the relation 


sin y 0 

♦1 " COS ( Y q - Y> *2 + *10 

where <j>io is a constant determined by initial conditions. 
Let the constant parameter w be defined as 


(26) 


s1n v 0 

W COS ( Yq " Y> 

Combining Eqs. (14), (25), and (27) leads to the relation 


(27) 




w 


(28) 


CIRCULAR CUTTER SURFACE 

In Fig. 9 the rotation of the cutter with a straight blade decribes a con- 
ical surface of revolution with vertex angle (it - 2yq) • The mean radius of 
the head cutter measured in the plane Z c = 0 is r c . The apex of the cone is 
at V with coordinates (0,0 , zq) in S c . The coordinates (x c , y c , z c ) of an 
arbitrary point C on the surface of revolution can then be expressed by the 
surface coordinates z c and a as 


x c = (z 0 ' z c )cot *0 cos 06 
y c = <z 0 " z c )cot *0 sin a 


(29) 
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Equation (29) may also be expressed In the form 

f(x c ,y c ,2 c > - tan 2 t Q (x 2 ♦ y 2 ) - (z Q - z c > 2 - 0 (30) 

The cutter surface may be expressed In the gear blank system S? by 
substituting from Eqs. (15), (20), and (22) Into Eq. (17). This leads to 

?c ■ t« c ,]tl« ) ,]tR 9|( HR p2 )? 2 * tR c| H V T sp ♦ tR c ,]T, c - to?* * T 

(31) 

where CE] and T are defined as 

(E) - tR cl HR )9 ]ER gp ]tR p2 3 (32) 


T -«cl”V T # P tt ''c1 ,T lc <33> 

By substituting from Eqs. (16), (18), (19), (23), and (24) Into Eqs. (32) and 
(33), the elements of CEJ and T, ejj and tfd.j, - 1,3), are found to be 


and 


e 12 

*13 

*21 

*22 

*23 

*31 

*32 

*33 


cos 4 2 cos (6 ♦ 4 ]) ♦ sin y sin 4 2 s 1 n(© ♦ ♦,) 
-sin + 2 cos (0 ♦ 4 j) ♦ sin y cos 4 2 s 1 n(© ♦ ♦ 1 ) 
-cos y s 1 n (6 ♦ 4 ^ ) 

-cos 4 2 s 1 n(© ♦ 4 ^) ♦ sin y sin 42 cos(© ♦ 4 ^> 

sin 4 2 s 1 n (6 ♦ 4 ^) ♦ sin y cos 4 2 cos (6 ♦ 4 1 > 

-cos y cos(© ♦ 4 j> 

cos y sin 4 2 

cos y cos 4 2 

sin y 



(34) 


» H cos © ♦ V sin © 
t 2 » -H sin © ♦ V cos © 
t, • -h 


(35) 
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Hence, the cutter surface expressed in the S 2 coordinate system may be 
obtained by substituting from Eqs. (31), (34), and (35) into Eq. (29). That 
is, 


(e ll 

tan 

4» 0 + 

e 31 cos a)x 2 + 

(e^tan 

*0 + 

e 32 cos a)y 2 + (e 13 tan 

^0 

+ e 33 COS 

a) Z 2 







= (Zq - t 3 )COS a - 

*1 

tan 


and 










(e 21 

tan 

+0 + 

e 31 si n a)x 2 + 

(e 12 tan 

*0 + 

e 32 sin a)y 2 + (e 13 tan 

*0 

+ e 33 sin 

a) Z 2 







= (z Q - t 3 ) s i n a - t 2 

tan 

^0 

( 36 ) 


SPIRAL BEVEL GEAR TOOTH SURFACE EQUATION 

Equations (36) represent the generating surface seen by the gear blank. 
In determining the envelope of the cutter surface on the gear blank, it is 

useful to evaluate the partial derivative of ? c in Eq. (31) with respect to 

the parameter <j> 2 - That is, 



3E 

34> 2 


*2 + 



3T 

94> 2 


(37) 


The second term is zero since r 2 is fixed in S 2 - From Eq. (35) the last 
term is also zero because the cutting speed is independent of the gear blank 
rotation. Let the component of the first term be djj. Then d}j can be 
determined from Eq. (34), and Eq. (37) may be rewritten in the form 

/ 


8r c 

^ - Cd, Jr 


3<J> 


i j 2 


(ij = 1,3) 


(38) 


where 

d 1 1 = -sin <|> 2 cos(0 + <t> ^ ) ( 1 - w sin y) + (sin y - w)cos <j> 2 
d 12 = -cos 4>2 cos(0 + (j^Kl - w sin y) - (sin y - w) sin <J> 
d 1 3 = -w cos y cos(0 + 4> 1 ) 

d 2 ^ = sin <j >2 sin(0 + <J^)(1 - w sin y) + (sin y - w) cos $2 
d 22 = cos 4> 2 s i n (0 + <^>(1 - w sin y) - (sin y - w)sin 4> 2 

^23 = ~ w cos Y + 4*1 ^ 

d 31 = cos y cos <j> 2 

d 32 = -cos y sin <j> 2 


s 1 n (0 + 4 ^ ) 

2 sin (0 + <j> 1 ) 


cos (0 + 4> 1 
cos (0 + ) 


J 


( 39 ) 


12 



Equation (38) represents the derivative of the equation giving transforma- 
tion from the cutter system to the gear blank system. It is useful for deriv- 
ing the constraint equation of the cutter's motion. The derivative of the 
cutter conical surface with respect to 4 > 2 , the rotation angle of the gear 
blank, is (from Eq. (30)) 


4 f( VW - 2 Un2 fo 


X c 9<j>, + y c 


3z 


+ 2(z rt - z .) 


0 ‘c' 9<J> 


(40) 


If z c is not equal to zq, we may substitute from Eqs. (29) and (38) 
i nto Eq. (40) . This leads to 

[tan ^(d^ cos a + d 21 sin a) + 1 + [tan Yq^i cos 06 + ^22 s ^ n a) 

+ d 32 ^2 + ^ an +0^13 C0S a + ^23 S ^ n a) + ^ 33 ^ Z 2 = 0 (41) 


Equation (41) is the constraint equation of the cutter motion on the gear 
blank. Combining Eqs. (36) and (41) forms a set of simultaneous equations 
representing the envelope of the cutter relative to the gear blank. The solu- 
tion of the simultaneous equations describes the tooth surface impression cre- 
ated by the cutter. Observe that since the coefficients are functions of (a, 
4 > 1 ) , the solution has the parametric form 


x 2 = x 2 (a,<J>i ) 

y 2 = y 2 (a, V 

z 2 = Z 2 (ot , 


(42) 


where a and <|>i are the surface coordinates 


DISCUSSION 

Equations (11) to (13) and Eqs. (36), (41), and (42) determine the cutter 
envelopes and hence the gear tooth surface for straight and spiral bevel 
gears, respectively. They form the basis for a numerical and computer graphic 
representation of the tooth surface. In this context these equations repre- 
sent an extension of the procedure described earlier (Chang, Huston, and Coy, 
1984) from spur gear to bevel gear. The difference here, however, is that the 
equations are more detailed and extensive. Hence, numerical and computer 
graphic analyses are needed. 

In a general sense the method simulates the kinematic relation of a 
cone-shaped gear blank and the cutter. With the rotation angle of the gear 
blank being the parameter, the envelope of the cutter profile on the gear 
blank describes the gear tooth profile. The straight and spiral bevel gear 
tooth surfaces are the forms given in Eqs. (13) and (42). The machine set- 
tings, the cutter radius, and the mean cone distance are the parameters that 
can be varied in a design optimization analysis. An accurate finite element 
mesh can also be obtained by discretizing the tooth surface. 
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In summary, it is believed by the authors that the analysis presented 
herein can form a basis for numerical design studies as well as for stress and 
deformation studies. Finally, the method may be readily extended to the study 
of nonintersecting shaft (hypoid) gears. 
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